Coevolution of Glauber-like Ising dynamics and topology 
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We study the coevolution of a generalized Glauber dynamics for Ising spins, with tunable thresh- 
old, and of the graph topology where the dynamics takes place. This simple coevolution dynamics 
generates a rich phase diagram in the space of the two parameters of the model, the threshold and 
the rewiring probability. The diagram displays phase transitions of different types: spin ordering, 
percolation, connectedness. At variance with traditional coevolution models, in which all spins of 
each connected component of the graph have equal value in the stationary state, we find that, for 
suitable choices of the parameters, the system may converge to a state in which spins of opposite 
sign coexist in the same component, organized in compact clusters of like-signed spins. Mean field 
calculations enable one to estimate some features of the phase diagram. 
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In recent times there has been an increasing attention, 
by the statistical physics community, towards applica- 
tions to social systems and relative phenomena The 
goal is the description and possibly the prediction of col- 
lective features of processes involving large numbers of 
individuals without detailed information on the charac- 
teristics of the single elements, much like it happens in 
the physics of phase transitions Q. Many simple mod- 
els have been devised, inspired by intuitive ideas on how 
social interactions between individuals take place. Such 
models are often variations of known models of statisti- 
cal physics, or entirely new and interesting types of dy- 
namics. The main ingredients are a graph, represent- 
ing the social network of interactions (acquaintances) 
between individuals, and a set of local rules, indicat- 
ing how the state of an agent is affected by (or affects) 
the state of its neighbors. The graph may be a lat- 
tice or have a more complex topology, reflecting prop- 
erties observed in real social networks [!, 0, [Hj]. Usu- 
ally one studies the model dynamics on a given graph 
topology, which remains frozen during the whole evo- 
lution of the process. However, in real social phenom- 
ena the dynamics of states is often coupled to the trans- 
formation of the social network where the process takes 
place, as the network evolves as well, and the time scales 
of the two evolutions may be comparable. So, a realis- 
tic description of social processes must consider the co- 
evolution of state dynamics and network topology. In 
the lastyears several coevolution models have been pro- 
posed 0, 0, H I, EE El EI EI EI EE EE Ezt The 



interaction rules of such models combine both changes in 
the states of the agents and in the link structure of the 
underlying graph. Frozen states of the dynamics are usu- 
ally characterized by a network composed of one or more 
connected components with all agents in each component 
being in the same state. Indeed, the dynamics of states 
in each component becomes independent of the dynam- 



ics ruling the states of the other components and, while 
the agents of each component converge to the same state, 
such a state will usually be different from one component 
to another. In several models both scenarios, i. e. one 
component with all agents in the same state and two or 
more separate components each in a different state, can 
be reached by suitable choices of the parameter weigh- 
ing the relative importance of the dynamics of the states 
versus that of the graph topology |& i, El El, EI El- 
Such a scenario is quite simple but it is not very realis- 
tic. For example, these models cannot describe the situa- 
tion in which different groups of people sharing the same 
state (domains) coexist in the same component, some- 
thing which is likely to happen in society. In this letter, 
we present the first model that accounts for this situation 
as well. Our model is based on a simple Glauber-type dy- 
namics for Ising spins [lij . It can also be seen as a sort 
of threshold model [IE] where disorder is in the topology 
and not in the thresholds. We show that, in spite of its 
simplicity, the model has a very rich behavior, with sev- 
eral phases, separated by transitions involving both the 
spin states and the graph topology. 

The starting point is a random graph a la Erdos- 
Renyi with N nodes and M links, with M = 
(k) being the average degree of the graph. We stress that 
the main results do not depend on the initial network 
topology because the rewiring dynamics leads inevitably 
to a random network with a Poisson degree distribution. 
Agents lie on the nodes of the graph, and are endowed 
with binary states (spins) a = ±1, which are initially 
assigned at random with equal probability 1/2. The dy- 
namics of the model is defined by iterating the following 
update rule: 

1. A node i is selected at random: we indicate with ki 
the number of its neighbors and with li the number 
of neighbors in the same state. 
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2. If li/ki > s, where < s < 1, the node is stable 
and nothing happens; otherwise a neighbor j with 
(jj 7^ Gi is randomly chosen and 

• with probability 0, i cuts its link to j and 

attaches it to a randomly chosen node / such 
that a\ = Gi and / is not already connected to 
i [IK (rewiring); 

• with probability 1 

flip). 



i adopts j's state (spin 



The model has two relevant parameters: the threshold s 
and the probability of rewiring (j). The threshold sets the 
minimum fraction of neighbors in the same state that 
a node must have to be stable. In this respect it is a 
measure of the sensitivity of agents against the social 
pressure exterted by neighbors with opposite state. If s 
is very small virtually all nodes are stable, i. e. they do 
not flip their spin nor rewire their connections. On the 
contrary, if s is close to 1 only nodes fully surrounded 
by nodes in the same state are stable. When s = 1/2, 
the spin dynamics is essentially the Glauber dynamics of 
Ising spins at zero temperature. For s < 1/2 the dynam- 
ics is rather uninteresting, as already in the initial condi- 
tion typically nodes have at least half of the neighbors in 
the same state and hence they are stable. Therefore we 
focus on the range 1/2 < s < 1. In general, the presence 
of a threshold allows for the existence of unsatisfied links 
(i.e. links joining agents with different states) in stable 
states of the system, at variance with standard revo- 
lution models. The rewiring probability <p is a measure 
of the relative importance of the rate at which the net- 
work evolves with respect to the rate at which the state 
of the node changes. The extreme values correspond to 
pure spin dynamics on a fixed network topology (0 = 0), 
and to pure network evolution, with no spin dynamics 

(0 = 1). 

The phase diagram of the model has a remarkably 
rich structure. To study its features, we monitor the 
behavior of some standard observables, the magnetiza- 
tion m = ^ZiCTi/N and the density of unsatisfied links 
p = Y^i<j — 5((Ji, Gj)]/M, where is the element 
of the adjacency matrix of the graph (A^ = 1, if i and j 
are neighbors, otherwise = 0), and 5 is Kronecker's 
delta function. Moreover we consider the convergence 
time t Cl defined as the time needed to reach a frozen con- 
figuration, at which the dynamics stops (i.e. k/ki > s 
for any i). In Fig. [I] we report in the plane (s, (j)) the 
numerical value of t Cl as well as the value of |ra| and p 
and the number of connected components of the graph, 
after a very long run, sufficient to reach the stationary 
state. From the behavior of the convergence time, it is 
clear that the parameter space is divided into two re- 
gions. In the upper left zone (denoted as S — 0) there 
is a phase with ongoing dynamic activity. In this region 
the convergence time t c diverges exponentially with iV, 



so that it is effectively infinite for systems of any rea- 
sonable size. However, the system reaches a stationary 
state, with constant value of the observables. Elsewhere, 
instead, the dynamics leads in a finite time to an ab- 
sorbing frozen state, with no dynamics. The two regions 
are separated by an absorbing-state phase transition [2l[ . 
The values of |ra| and p indicate that the active phase 
is disordered: the average magnetization remains zero 
and the density of unsatisfied links remains high. In this 
phase, due to the high value of s, sites are rarely stable 
and they keep rewiring, looking for similar partners, as 
m Refs. [11] and 0. The structure of the absorbing 



phase is much richer, as one can identify several distinct 
subphases, with various types of internal organization. 
For small values of (j) and s, there is a phase (A-I) where 
| ra | = 1 and p = 0. Here the relatively slow rewiring 
process allows spin ordering to be completed while the 
topology remains globally connected. The opposite oc- 
curs in the upper right corner of the parameter space 
(phase A-II). For large s and 0, the value of p is the 
same of the A-I zone, but in the A-II zone the number 
of connected components is 2: the system splits into two 
topologically separated sets of similar size, one of them 
fully ordered with g = 1 and the other with g = — 1. 
Phases A-I and A-II correspond to those found in the 
voter model . To understand the organization of the 
system in the rest of the plane, we study the presence and 
extent of homogeneous domains, intended as subsets of 
nodes of the network with two properties: 1) all nodes of 
the subset are in the same state and 2) any pair of nodes 
of the subset can be joined with a path within the sub- 
set. We then measure two new observables: the number 
of homogeneous domains in the network and the relative 
size of the largest of them. In Fig. [2] we plot these two 
observables as a function of s for a fixed large value of 
the rewiring probability <p. It is possible to identify two 
new phases, delimited by two threshold values s p and s q 
(indicated by the two vertical dotted lines). 

For s < s p there are just microscopic domains and 
their number is proportional to the number of nodes in 
the network. In this phase (A — 0), that spans the whole 
range of for s < 1/2, ra = and p ~ 1/2. Stability 
is rather easy to reach for all nodes, after few spin flips 
or link rewirings. We stress that A — and 5 — are 
different: in the case of A — an absorbing state is al- 
ways reached, while for S — the system reaches only a 
dynamic stationary state. At s = s p a percolation tran- 
sition [22|] takes place: for s p < s < s q the dynamics 
lasts long enough to allow for the formation of macro- 
scopic domains (typically two) that grow bigger as s is 
increased up to the point where they occupy the whole 
system (phase A — D). Finally, for still higher values 
of the threshold s > s q the macroscopic domains be- 
come topologically disconnected from each other and co- 
incide with the two connected components of the network 
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FIG. 1: (color online) Dependence on the threshold s and the rewiring probability <j) of several variables: the convergence time 
t c (top left), the absolute value of the magnetization (top right), the density of unsatisfied links (bottom left) and the number 
of connected components (bottom right). Data are obtained from numerical simulations on a graph with N = 50000, (k) = 10. 
The darker color in the top left panel means that the convergence time is larger than 150 times the number of nodes. Data in 
the three other panels are computed after 1.5 x 10 6 iterations. 



(phase A — II). Based on this evidence, we schematically 
represent in Fig. [3] the phase-diagram of the model. Our 
simulations, performed up to size TV = 50000, seem to 
indicate that the parameter space is divided into genuine 
phases separated by well-defined transition lines. How- 
ever, a detailed investigation of the nature of all of them 
(and of the associated critical behavior) is numerically 
very demanding and goes beyond the scope of the present 
paper. 



Some of the features of the phase-diagram are recov- 
ered (see Fig. [3]) via a mean field (MF) approach, similar 
to the one in Ref. [13]. At each step, a node with k links 
is randomly chosen. Since the rewiring dynamics leads to 
a network with a Poisson degree distribution, the proba- 
bility Pk(t) to extract a node with k links is supposed to 
be a Poisson distribution with mean (k) at each step. De- 
note with j the number of unsatisfied links and with Vkj 
the probability that the chosen node is not stable and 
hence must be updated. With probability (j)Vk,j a ran- 



dom unsatisfied link is rewired. In this case, the density 
of unsatisfied links changes by Ap = —j^- On the con- 
trary, with probability (1 — 0) Vkj the state of the node 
is flipped and the density of unsatisfied links changes by 
the quantity Ap = ^j^p- • Using these expressions, it is 
possible to write the time evolution master equation for 
a generic update rule 



dp 



W = V — V B ip) 
dt k 'i 



/N 



Vk, : 



3=0 



(1-0) 



2 (k - 2j) 
(k)N 



(k)N m 
(1) 



where 1/N is the temporal interval between successive 
steps and B^- is the probability to find a node with j 
unsatisfied links. In a MF spirit, all nodes of the network 
can be considered equivalent and the probability to have 
an unsatisfied link can depend only on the global observ- 
able p. Thus, the probability to have an unsatisfied link 
is taken independent for each node and is well approxi- 
mated by a binomial distribution B^ = p 3 ; (1 — p) k ~i . 
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FIG. 2: (color online) Number of domains and density of the 
largest domain as a function of s, obtained from simulations 
on a graph with (k) = 10 and (f> = 0.95. 
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the right-hand side of Eq. [T] has the simple expression 
F^(p) = 2(l-0)(l-2p) 

-e-< fc >^2(l-0)(l-p)- ( 2 ) 
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m— For p s sufficiently small, the stationary solution of Eq. (pQ) 



2[(fe)(l-0)-l] 



For 



< 



has then the form p s - (fc)[2+(fc)(1 _ 0) _ 30 ] • 
c (s = 1) = there is an active stationary state with 

p s > 0. For larger values of qb the density of unsatisfied 
links is zero, corresponding to an absorbing phase. The 
transition is predicted to be continuous. For (k) = 10 
the critical value is (j) c (s = 1) = 0.9, in agreement with 
numerical simulations. It is possible to determine numer- 
ically the transition line (j) c (s) for any value of s. The 
resulting curve is reported in Fig. [3j 

The simple model we have proposed offers a surpris- 
ingly rich variety of possible scenarios, by varying the 
two parameters s and (j). In particular, phase boundaries 
correspond to magnetization, connectedness and/or per- 
colation transitions. The most striking feature, absent 
in all other models of coevolution, is the existence of a 
phase where stable homogeneous domains coexist in the 
system, even if the latter is not split into components. 
This feature is due to the presence of the threshold s: 
models characterized by a threshold are likely to display 
this type of behavior and represent a promising option 
for a realistic description of social phenomena. We stress 
however that the goal of this paper was not a description 
of a specific real world phenomenon, but rather the in- 
vestigation of what are the possible qualitative outcomes 
when threshold dynamics and rewiring operate simulta- 
neously. 



FIG. 3: (color online) Schematic representation of the phase 
diagram obtained from numerical simulations on a graph with 



N = 50000, (k) 



10, after 1.5 x 10 iterations. We can 



identify two types of phases: the dynamic stationary phase 
(denoted by S — 0) and four absorbing phases (denoted by 
A), divided by the black solid line. Each absorbing phase 
is characterized by a different organization of the network 
in domains (see text for details). The (red) dot-dashed line 
corresponds to the numerical solution of mean-field equations. 



In the case of the voter model the update proba- 
bility is simply j/k, and it is possible to analytically 
solve it [13[. In our model, the update probability is 
*Pkj — Q (j/k + 5 — 1) where (x) is the Heaviside step 
function. Due to the nonlinearity of Eq. ([1]) an analytical 
expression can be found only for s = 1. In such a case 
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